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Abstract 

A  wavelet-collocation  scheme,  constructed  from  the  discrete  singular  convolution  (DSC) 
is  introduced  for  computational  electromagnetics.  The  basic  philosophy  behind  the  DSC 
algorithm  is  studied.  Four  test  problems,  including  waveguide  analysis  in  both  regular  and 
irregular  domains,  electromagnetic  wave  scattering,  electrostatic  field  estimation  via  poten¬ 
tial  functions,  and  electromagnetic  wave  propagation  in  three  spatial  dimensions,  are  em¬ 
ployed  to  illustrate  the  usefulness,  test  the  accuracy  and  explore  the  limitation  of  the  wavelet 
algorithm.  Computational  accuracy  of  the  wavelet  algorithm  is  shown  to  be  controllable  and 
is  compared  with  a  standard  approach.  Numerical  experiments  indicate  that  the  proposed 
wavelet  algorithm  is  a  promising  approach  for  solving  problems  in  electromagnetics. 

Keywords:  Maxwell’s  equations,  waveguide,  wave  scattering,  potential  function,  wave  propaga¬ 
tion,  wavelet-collocation,  discrete  singular  convolution. 

1  Introduction 

Recently,  computational  electromagnetics  (GEM)  has  emerged  as  a  distinct  scientific  discipline  for 
the  study  and  understanding  of  a  wide  variety  of  electrical  and  electronic  engineering  problems  [1, 
2,  3,  4,  5,  6,  7,  8].  As  a  natural  extension  to  the  analytical  approach,  the  GEM  is  based  on 
numerically  solving  the  governing  equations  either  in  the  partial  differential  form  or  in  the  integral 
equation  form.  The  complexities  of  physics  and  of  the  field  geometry  are  no  longer  the  limiting 
factors  to  GEM  as  they  are  to  the  analytical  approach.  With  the  advent  of  high-performance 
digital  computers,  GEM  is  emerging  as  a  powerful  approach  for  solving  practical  problems  in 
electromagnetics  [1,  2],  Asa  computational  discipline,  the  success  of  the  GEM  is  vitally  dependent 
on  the  virtues  of  computational  algorithms,  such  as  numerical  accuracy,  stability  and  efficiency. 
These  in  turn  depend  on  grid  methods  and  numerical  schemes  for  the  solution  of  the  Maxwell’s 
equations. 

Typically,  grid  methods  used  in  GEM  are  either  local,  such  as  finite  difference  (FD)  methods 
[2],  or  global,  such  as  Fourier  pseudospectral  methods  [3].  Global  methods  are  highly  accurate 
but  are  cumbersome  to  implement  in  complex  geometries  and  non-conventional  boundary  condi¬ 
tions.  In  contrast,  local  methods  are  easy  to  implement  for  complex  geometries  and  discontinuous 
boundary  conditions.  However,  the  accuracy  of  local  methods  is  relatively  low[4].  There  exists 
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many  problems  in  CEM  which  require  both  high  computational  accuracy  and  high  numerical  flex¬ 
ibility  in  handling  complex  geometries.  These  problems  are  characterized  by  a  geometry  which 
has  a  large  domain  size,  i.e.,  the  dimensions  of  the  scatterer  greatly  exceed  the  wavelength  of  the 
incident  wave.  A  typical  example  is  the  radar  cross-section  analysis  of  an  entire  airplane  with 
an  incident  electromagnetic  wave  having  a  frequency  of  the  order  of  ten  GHz.  To  deal  with  such 
problems,  it  is  desired  to  have  a  computational  method  that  has  both  global  methods’  accuracy 
and  local  methods’  flexibility. 

Wavelet-based  algorithms  are  expected  to  fill  this  gap,  and  wavelet  analysis  has  been  applied 
to  CEM  [7,  8].  More  recently,  the  discrete  singular  convolution  (DSC)  algorithm  was  proposed  as 
a  potential  wavelet  approach  for  the  computer  realization  of  singular  convolutions [9].  The  math¬ 
ematical  foundation  of  the  algorithm  is  the  theory  of  distributions  and  the  theory  of  wavelets 
[10,  11].  By  appropriately  selecting  parameters  of  a  DSC-wavelet  kernel,  the  DSC-wavelet  ap¬ 
proach  exhibits  controllable  accuracy  for  integration  and  shows  great  flexibility  in  handling  com¬ 
plex  geometries  and  boundary  conditions.  It  was  demonstrated  that  in  the  framework  of  the  DSC 
algorithm,  different  numerical  schemes,  such  as  global,  local,  Galerkin,  collocation,  and  FD,  can 
be  deduced  from  a  single  starting  point  [12], 

The  purpose  of  this  paper  is  to  introduce  the  DSC-wavelet  algorithm  for  CEM.  The  potential 
of  the  wavelet  algorithm  for  CEM  is  explored  through  applications  to  waveguide  analysis,  wave 
scattering,  electrostatics  analysis,  and  wave  propagation. 

2  The  Discrete  Singular  Convolution  and  wavelets 

In  the  DSC-wavelet  algorithm,  it  is  convenient  to  formulate  the  numerical  solution  in  terms  of 
the  collocation  scheme  for  solving  PDEs  [12],  For  the  collocation  scheme,  an  approximation  to 
the  derivatives  is  required  given  that  function  values  {f(xm)}  at  a  set  of  discrete  points  {xm}  are 
known.  In  the  DSC  algorithm,  we  approximate  the  function  and  its  nth  order  derivative  at  point 
x  by  a  discrete  convolution 

M 

f(n\x)  «  5$r(x~xk)f(xk)  (n  =  o,  1,2,  •  ■  ■),  (1) 

k=—M 

where  2 M  +  1  is  the  computational  bandwidth  and  {x^  are  centred  around  x.  Here  8ajCr(x  —  Xk) 
is  a  collective  symbol  for  one  of  the  (regularized)  DSC  kernels,  or  in  general,  any  delta  sequence 
kernel  which  provides  an  approximation  to  the  delta  distribution  5  [9].  For  a  given  5a>a(x  —  Xk ), 
the  higher  order  derivative  terms  can  be  given  by  means  of  analytical  differentiation  S^A(x  —  Xk)  = 
[{£)n^Ax  ~~  xk)\-  Numerical  solution  of  differential  equations  then  can  be  easily  implemented 
in  a  collocation  scheme  by  using  Eq.  (1).  Moreover,  the  wavelet  collocation  has  been  shown  to 
be  equivalent  to  a  wavelet  Galerkin  [12], 

In  the  DSC  algorithm,  delta  sequence  kernels  can  be  constructed  as  either  positive  type  or 
Dirichlet  type  [11].  The  Dirichlet  type  kernels  are  well-suited  to  be  utilized  in  a  wavelet-collocation 
approach,  while  the  positive  type  kernels  can  be  used  in  a  wavelet- Galerkin  scheme.  A  family 
wavelet  delta  kernels  including  both  types  have  been  constructed  [9].  For  simplicity,  only  a  few 
relevant  delta  sequence  kernels  of  Dirichlet  type  are  considered  here,  i.e.,  Dirichlet  kernel,  Lagrange 
kernel,  and  Shannon  kernel.  The  Shannon  kernel  is  the  well  known  scaling  function  for  Shannon 
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wavelets,  which  play  a  crucial  role  in  information  theory  and  the  theory  of  signal  processing. 
However  the  usefulness  of  Shannon  delta  kernel  for  solving  PDE  is  limited  by  the  fact  that  it  is 
non-local. 

In  the  DSC-wavelet  algorithm,  a  regularization  procedure  is  typically  applied  to  improve  the 
smoothness,  regularity  and  localization  of  a  temper  distribution  by  using  the  Gaussian  [9].  In 
this  manner,  we  construct  the  following  three  DSC-wavelet  kernels,  i.e.,  a  regularized  Shannon’s 
kernel  (RSK),  a  regularized  Dirichlet  kernel  (RDK),  and  a  regularized  Lagrange  kernel  (RLK)  [9], 
where  A  is  the  grid  spacing, 


Sa, <t(x  -  xk) 


sin  f  (x  -  xk) 

Ux~xk) 


exp 


(x  -  xk )2 
2a2 


(2) 


SA,a(X  ~  Xk) 


sin  [f  (x  -  xk ) 

(2^  +  l)sin[iEff]  °XP 


(x  -  xk )2 
2<72 


Sa Ax  ~  xk) 


i=k+L 

n 

i=k-L,i^k 


Xk  Xi 


exp 


(x  -  Xk)2 
2a 2 


(3) 

(4) 


3  Numerical  experiments 

In  this  section,  we  explore  the  utility  of  the  wavelet-collocation  approach  for  solving  problems  in 
electromagnetics.  Four  different  GEM  problems  are  considered.  In  the  wavelet  algorithm,  the  half 
bandwidth,  M,  is  set  to  40  in  all  calculations  except  where  specified.  In  the  regularized  Dirichlet 
and  Lagrange  kernels,  (3)  and  (4),  L  is  set  to  50  in  calculations.  The  wavelet  parameter  r  =  a/ A 
is  optimized  for  each  M  and  reported.  The  standard  second-order  central  FD  method  is  also 
considered  for  a  comparison. 


3.1  Waveguide  analysis 

Eigenmode  ANALYSIS  We  first  consider  the  problem  of  finding  waveguide  eigen-modes.  Consider 
a  rectangular  waveguide,  with  propagation  in  the  ^-direction.  For  the  transverse  magnetic  (TM) 
mode,  the  Ez  field  component  can  be  expressed  as  Ez(x,  y ,  z}  t )  =  E(x,  where  E  satisfies 

the  Helmholtz  equation 


d2E 

dx2 


d2E 

dy2 


k2E  =  0, 


(5) 


and  vanishes  on  the  boundary  ( E  =  0)  for  the  TM  modes.  The  eigenvalue  k 2  determines  the 
phase  parameter  a  through  k2  =  u2ye  —  a2,  where  e  and  ji  are  dielectric  constant  and  magnetic 
permeability,  respectively. 

In  numerical  computation,  different  grid  size  N  is  considered  for  a  computational  domain, 
[0, 107r]  x  [0, 107r],  and  the  anti-symmetric  boundary  extension  scheme  is  employed  for  the  PEC 
boundaries  [9].  By  considering  both  the  wavelet  and  FD  methods,  the  numerical  results  are  listed 
in  Table  1.  Here  for  three  DSC-wavelet  kernels,  we  choose  M  =  N .  When  bandwidths  are 
M  =  36,24  and  12,  the  values  oir  —  a/A  are  chosen  as  4.2,  3.2  and  2.65  for  the  RSK  and  RDK 
and  2.95,  2.75  and  2.3  for  the  RLK.  It  is  obvious  in  Table  1  that  three  DSC-wavelet  schemes 
are  significantly  more  accurate  than  the  FD  method,  especially  when  representing  high-order 
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Table  1:  Absolute  errors  for  waveguide  eigenmode  analysis.  Here  1.7(-2)  denotes  1.7  x  10 


-2 


N 

Scheme 

Eigenmode 

No.  10 

No.  20 

No.  40 

No.  60 

No.  80 

No.  100 

FD 

1.7(-2) 

6.5(-2) 

2.5(4) 

3.9(-l) 

6.4(4) 

8.6(4) 

1  9 

RSK 

8.4(-8) 

5.2(-7) 

1,1(4) 

1.9(-1) 

3.2(4) 

4.2(4) 

JL  Z 

RDK 

8.8(-8) 

5.1(-7) 

1.1(4) 

1.9(-1) 

3.2(4) 

4.2(4) 

RLK 

2-l(-7) 

4.7  (-7) 

1.1(4) 

1.9(-1) 

3.2(4) 

4.2(4) 

FD 

4.0(-3) 

7.9(-3) 

2.9(-2) 

5.6(-2) 

1.5(-1) 

3.0(-l) 

0  A 

RSK 

1.8(-14) 

9.5(-15) 

8.7(-13) 

1.4(-8) 

8.0(-8) 

6.4(-9) 

RDK 

3.1(-15) 

1.6(-14) 

8.5(-13) 

1.3(-8) 

7.8(-8) 

6.3(-9) 

RLK 

9.7(-13) 

3.2(-12) 

5.1(-13) 

3.5(-10) 

6.3(-9) 

3.6(-10) 

FD 

1.7(-3) 

3.4(-3) 

1.3(-2) 

2.5(-2) 

6.7(-2) 

1.3(-1) 

RSK 

1.2(-14) 

4.3(-14) 

2.2(-16) 

7.4(-14) 

2.2(-14) 

3.7(-14) 

ou 

RDK 

2.1(-14) 

2.1(-14) 

3.9(-14) 

2.0(-14) 

4.2(-14) 

1.8(-15) 

RLK 

2.9(-14) 

1.2(-15) 

5.0(-15) 

3.8(-14) 

6.7(-15) 

3.3(-14) 

eigenmodes.  Furthermore,  we  note  that  the  RSK  and  RDK  have  almost  identical  performance  in 
all  computations.  Hence,  we  only  need  to  consider  one  of  them  in  the  rest  of  the  computation. 
The  performance  of  the  RLK  differs  slightly  from  that  of  the  RSK.  However,  such  differences  are 
extremely  small.  Thus,  one  can  argue  that  each  wavelet  kernel  performs  equally  well.  This  finding 
is  in  consistent  with  those  in  [9,  12], 

For  N  =  12,  there  are  only  a  maximum  of  2.4  grid  points  per  wavelength  (PPW)  in  each 
dimension.  Obviously,  if  this  is  the  case,  it  is  impossible  for  any  numerical  method  to  exactly  catch 
up  with  the  fast  oscillations  in  high-order  eigenmodes.  However,  under  such  a  low  grid  density,  the 
DSC-wavelet  results  are  still  much  better  than  those  of  the  FD.  All  three  wavelet  kernels  deliver 
excellent  results  when  the  mesh  is  refined  to  242  (about  4.8  PPW  in  each  dimension).  Moreover, 
machine  precision  is  attained  for  the  first  100  eigenmodes  when  the  mesh  is  further  refined  to 
362.  However,  the  accuracy  improvement  of  the  FD  scheme  is  negligible  and  the  final  results  for 
high-order  eigenmodes  remain  inaccurate. 

Field  distribution  In  practical  applications,  the  geometry  of  a  waveguide  can  be  very 
complex  so  as  to  achieve  certain  required  properties.  Therefore,  it  is  important  for  a  computational 
algorithm  to  be  flexible  enough  in  handling  complex  structures.  To  test  the  wavelet  algorithm 
in  this  regard,  we  consider  a  specially  designed  waveguide  with  its  geometry  having  a  T-shape. 
We  are  interested  in  the  field  distribution  of  TM  (and/or  TE)  modes  of  the  waveguide.  By  using 
the  RSK  with  r  =  4.8,  the  obtained  field  distributions  across  the  cross  section  for  eigenmodes  1, 
2,  3  and  4  are  shown  in  FIG.  1(a).  Due  to  complexity  in  geometry,  these  eigenmodes  are  quite 
localized  and  do  not  give  the  full  shape  of  waveguide  confining  geometry.  However,  as  plotted 
in  FIG.  1(a),  higher-order  eigenmodes,  17,  18,  19  and  20,  clearly  reflect  the  waveguide  geometry. 
The  confining  geometry  of  the  waveguide  can  be  seen  from  mode  patterns. 
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Figure  1:  Plots  of  eigenmodes  and  scattering  fields,  (a)  Eigenvectors  of  the  T-shape  waveguide 
obtained  by  using  the  RSK  and  a  mesh  size  of  50  x  50.  In  (a),  left:  eigenmodes  1,  2,  3  and  4;  right: 
eigenmodes  17,  18,  19  and  20.  (b)  Plots  of  the  real  part  of  the  scattering  wave  u.  Diamonds:  the 
exact  solution;  solid  lines:  the  DSC  results;  and  dashed  lines:  the  FD  results. 

3.2  Wave  scattering 

Classical  problems  of  wave  scattering  are  usually  modeled  through  an  exterior  problems  governed 
by  the  Helmholtz  equation  with  the  Sommerfeld  radiation  conditions  imposed  at  infinity  [5].  For 
the  numerical  solution  of  such  an  exterior  problem,  an  artificial  boundary  is  frequently  introduced 
[5].  In  ID,  the  resulting  model  problem  over  a  finite  region  can  be  given  as  [6] 

u"  +  k2n2(x)u  =  0,  in  [0,7r],  (6) 

with  'it(O)  =  1  and  u'( tt)  —  iku( tt)  =  0.  Here  k  is  the  wavenumber  and  n(x)  is  the  index  of 

refraction. 

By  choosing  k  —  25  and  n(x)  =  exp[— 10(x  —  l)2]  +  1  [6],  both  the  RSK  and  FD  schemes 
are  employed  to  solve  the  boundary  value  problem  (6).  In  the  RSK,  r  is  also  chosen  as  4.8.  By 
using  N  =  60,  240,  and  960,  the  real  parts  of  the  simulated  waves  are  plotted  in  FIG.  1(b).  In 

order  to  benchmark  the  numerical  results,  an  exact  solution  is  also  shown  in  FIG.  1(b)  which 

is  actually  an  accurate  approximation  obtained  by  using  the  FD  scheme  with  3840  points.  In 
terms  of  the  minimal  wavelength,  the  grid  density  of  this  problem  is  2.4,  9.6,  and  38.5  PPW 
for  N  =  60,  240,  and  960,  respectively.  Since  about  18  PPW  is  required  for  the  FD  scheme  for 
accurate  modeling  [4],  the  FD  results  are  graphically  accurate  only  when  N  =  960.  Amazingly, 
the  wavelet  approximation  is  indistinguishable  from  the  exact  solution  even  when  the  grid  density 
is  as  low  as  2.4  PPW.  Such  a  under  sampling  problem  remains  challenging  to  other  local  methods, 
while  the  mixed  boundary  condition  may  limit  the  applicability  of  global  methods. 

3.3  Electrostatics  analysis 

In  many  electromagnetic  applications,  potential  function  provides  a  useful  and  convenient  alter¬ 
native  representation  of  electromagnetic  fields.  For  electrostatic  problems,  formulation  with  the 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


Figure  2:  Contour  plots  of  the  electric  scale  potential  V.  (a)  Without  charge  density;  (b)  with  an 
area  of  charge  density. 


electric  scalar  potential  is  known  to  be  the  most  efficient  way  to  attain  the  solutions  [1].  By 
representing  the  electric  held  E  as  the  gradient  of  a  scalar  potential  E  =  —  VE,  on  attains  the 
electrostatic  equation  [1] 

VV  =  -P-i  (7) 

where  p  is  the  volume  electric  charge  density. 

To  demonstrate  the  wavelet  algorithm  for  electrostatic  calculations,  we  first  consider  a  simple 
problem  without  charge  density  (p  =  0).  The  geometry  of  the  problem  consists  of  a  conducting 
square  box  of  lm  x  lm  in  the  x  —  y  plane.  The  box  is  infinitely  long  in  the  ^-direction  and 
thus,  this  can  be  regarded  as  a  two-dimensional  problem.  The  potential  at  the  bottom  side  is 
considered  to  be  isolated  from  all  other  sides  with  values  V(x,  0)  =  10G  and  V(0,y)  =  V(l,y)  = 
V{x,l)  =  0.  By  choosing  32  points  (N  =  32)  in  each  dimension,  the  Laplace  equation  (7) 
for  scale  potential  function  can  be  directly  solved  by  using  the  DSC  algorithm  with  M  —  32 
and  r  =  3.8.  The  anti-symmetric  extension  is  used  for  all  boundaries.  The  results  are  plotted 
in  FIG.  2a.  We  next  introduce  charge  density  in  the  above  problem,  while  the  domain  and 
boundary  conditions  are  kept  the  same.  Inside  the  box,  an  area  of  charge  density  is  located  at 
x  :  0.41m  —  0.59m,  y  :  0.72 m  —  0.88m.  The  charge  density  is  uniformly  distributed  within  a 
material  of  relative  permittivity  er  =  25  and  has  a  value  of  p  =  1.0  x  10~'[C/m2].  Similar  to 
the  first  problem,  the  DSC-wavelet  approximations  of  this  inhomogeneous  problem  can  be  easily 
generated,  see  FIG.  2b. 


3.4  Electromagnetic  wave  propagation 


In  an  uncharged  homogeneous  media,  propagation  of  electromagnetic  waves  is  governed  by  the 
wave  equation 


d  2W 
d  t2 


—V2W  =  0, 

ep 


(8) 
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Table  2:  Maximum  errors  for  solving  the  3D  wave  equation. 


N 

Scheme 

M 

At 

t  =  4 

t  =  10 

t  =  16 

t  =  22 

FD 

— 

0.2 

S-l(-l) 

1.2(+0) 

1.7(+0) 

2.0(+0) 

6 

0.5 

l.l(-3) 

2.3(-3) 

3.5(-3) 

4.8(-3) 

RSK 

12 

0.2 

2.8(-5) 

6.1(-5) 

9.4(-5) 

1-3  (-4) 

24 

24 

0.05 

7.6(-8) 

1.4(-7) 

2.1(-7) 

2.7(-7) 

6 

0.5 

1.5(-3) 

3.3(-3) 

5.2(-3) 

7.0(-3) 

RLK 

12 

0.2 

3.1(-5) 

7-1  (-5) 

l.l(-4) 

1.5(-4) 

24 

0.05 

2.0(-7) 

4.9(-7) 

7.8(-7) 

l.l(-6) 

FD 

— 

0.2 

2'2(-l) 

5.6(-l) 

9.1(-1) 

1.2(+0) 

12 

0.2 

4.8(-5) 

l-2(-4) 

2.0(-4) 

2.6(-4) 

RSK 

24 

0.05 

2.1(-7) 

5.2(-7) 

8.2(-7) 

l.l(-6) 

36 

36 

0.001 

3.0(-13) 

3.5(-14) 

3.3(-12) 

4.1(-12) 

12 

0.2 

4.7(-5) 

1.2(-4) 

1.9  (-4) 

2.5(-4) 

RLK 

24 

0.05 

2.1(-7) 

5.2(-7) 

8.2(-7) 

l.l(-6) 

36 

0.001 

3.0(-13) 

2.7(-14) 

3.3(-12) 

4.1(-12) 

where  W (r,  t)  can  be  either  the  electric  held  or  the  magnetic  held.  To  illustrate  the  potential  of 
the  wavelet  algorithm  for  electromagnetic  wave  simulations,  we  consider  an  analytically  solvable 
problem 


W  (. x ,  y,  z,  t )  =  sin 


cixx  T  Oiyy  olzz  T 


^a2x  +  a2y  +  a2J(ey)t 


with  periodic  boundary  condition.  To  further  simplify  notations,  we  set  ax  =  ay  =  az  =  ey  =  1 
and  focus  only  on  a  cubic  domain  of  size  (107r)3. 

Two  wavelet  schemes:  the  RSK  and  RLK,  and  the  FD  method  are  considered  for  this  problem. 
To  test  the  accuracy  and  rate  of  convergence  of  the  wavelet  schemes,  a  variety  of  grid  points  N 
(N  =  Nx  =  Ny  =  Nz)  and  computational  bandwidth  M  (M  =  Mx  =  My  =  Mz)  are  employed. 
When  M  =36,  24,  12,  and  6,  the  values  of  r  are  4.2,  3.2,  2.65  and  2.0  for  the  RSK  and  2.95  2.75, 
2.3  and  1.85  for  the  RLK.  The  fourth  order  Runge-Kutta  scheme  is  used  in  all  schemes  for  time 
discretization.  Sufficiently  small  time  increments  (At)  are  used  so  that  the  major  errors  are  due  to 
the  spatial  discretization.  The  L ^  errors  at  a  number  of  time  units  are  listed  in  TABLE  2.  When 
N  =  24,  the  number  of  PPW  is  4.8  in  each  dimension.  This  is  a  typical  case  of  under  sampling 
and  it  is  very  difficult  to  achieve  high  computational  accuracy  by  any  means.  Since  usually  the 
Yee  algorithm  uses  about  18  PPW  [4] ,  the  FD  results  are  incorrect  when  N  =  24.  ffowever,  errors 
in  the  wavelet  results  for  this  difficult  problem  are  all  smaller  than  0.5%.  Moreover,  the  wavelet 
accuracy  significantly  improves  as  the  M  value  is  increased  from  6  to  12  and  24.  The  last  case, 
M  =  N  =  24,  corresponds  to  a  global  collocation  method  and  its  results  are  significantly  better. 
In  the  case  of  N  =  36,  it  is  on  an  average  about  7.2  PPW  in  each  dimension.  The  accuracy  of 
the  FD  method  is  still  unsatisfactory,  whereas  the  wavelet  solutions  are  more  accurate  than  those 
of  N  =  24.  Similarly,  the  accuracy  of  the  wavelet  results  dramatically  increases  as  M  increases. 
Furthermore,  it  is  interesting  to  note  that  the  wavelet  algorithm  attains  extremely  high  accuracy 
at  the  global  limit  of  M  =  N  =  36. 
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4  Conclusion 

In  conclusion,  this  paper  introduces  a  wavelet-collocation  scheme  based  on  the  discrete  singular 
convolution  (DSC)  algorithm  for  computational  electromagnetics  (CEM).  The  computational  phi¬ 
losophy  of  the  DSC- wavelet  algorithm  is  studied.  Three  typical  wavelet  kernels  are  utilized  to 
illustrate  the  usefulness  of  the  proposed  algorithm.  The  potential  of  the  wavelet  algorithm  for 
CEM  is  explored  by  studying  problems  of  waveguide  eigenmode  analysis,  wave  scattering  simula¬ 
tion,  electrostatics  analysis,  and  electromagnetic  wave  propagations.  The  standard  central  finite 
difference  (FD)  scheme  is  also  employed  for  a  comparison. 

It  is  found  that  the  wavelet  errors  are  usually  smaller  than  those  of  the  FD  method  by  several 
or  over  ten  orders  of  magnitude.  Especially,  the  wavelet  accuracy  is  found  to  be  superior  for 
numerically  challenging  problems  of  under  sampling.  The  wavelet  algorithm  is  applied  to  handle 
complex  geometries,  such  as  a  T-shape  waveguide  and  a  non-uniform  electrostatics  domain  with 
charge  density.  The  non-reflecting  boundary  condition  raised  in  wave  scattering  has  also  been 
successfully  treated  in  the  wavelet  algorithm.  In  summary,  present  work  indicates  that  the  wavelet 
algorithm  is  a  promising  and  potential  approach  for  computational  electromagnetics. 
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